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O I Current cosmological data drop an interesting hint about the existence of sub- 

"^ ' eV sterile neutrinos, which should be a part of the cosmic neutrino background 

Ci . (Ci^B). We point out that such light sterile neutrinos may leave a distinct imprint 

on the electron energy spectrum in the capture of relic electron neutrinos by means 
CN ■ of radioactive beta-decaying nuclei. We examine possible signals of sterile neutrinos 

relative to active neutrinos, characterized by their masses and sensitive to their 
number densities, in the reaction fg + ^H — )• ^He + e~ against the corresponding 
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On ■ tritium beta decay. We stress that this kind of direct laboratory detection of the 

^^ . Cz^B and its sterile component might not be hopeless in the long term. 
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1 1 1 As fairly stable and weakly interacting particles, relic neutrinos of the Big Bang 
must survive today and form a cosmic background similar to the cosmic microwave back- 
ground (CMB) radiation. This cosmic neutrino background (Ci^B) played an important 
role in the evolution of the Universe, and its existence has been indirectly "seen" from 
current cosmological data on the Big Bang nucleosynthesis (BBN), CMB anisotropies and 
large-scale structures [1]. How to directly detect the Cz/B in a laboratory experiment is a 
great challenge to the present experimental techniques, simply because today's tempera- 
ture of the Ci/B is extremely low (T^, ^ 1.945 K) and thus the average three-momentum of 
each relic neutrino is very small {{Py) = 3T^ f» 5 x 10~^ eV). Among several possibilities 
for the direct Cz/B detection [2], the most promising one seems to be the relic neutrino 
capture experiment by means of radioactive beta-decaying nuclei [3] — [8]. The point is 
that a generic neutrino capture reaction 1^^ + N ^ N' + e~ will take place with no thresh- 
old on the incident neutrino energy, provided A^ can naturally undergo the beta decay 
N ^ N' + e~ +17^ with an energy release Qa = rrij^ — m^, — m^ in the limit of vanishing 
neutrino masses (i.e., m,^ — )■ for z = 1, 2, ■ ■ ■). The signal of this neutrino capture process 
is measured by the monoenergetic electron's kinetic energy Q^ + E^^ > Qjs + rn^ for each 
neutrino mass eigenstate u^, as compared with the non-monoenergetic electron's endpoint 
energy Q^ — m^ for each u^ in the corresponding beta decay. So there is a gap equal to or 
larger than 2m- between the kinetic energies of the detected electrons in ly^ + N ^ N' + e~ 
(signal) and N ^ N' + e" + i7^ (background). A measurement of this gap will directly 
probe relic neutrinos and determine or constrain their masses. 

An immediate question is whether the Cz/B consists of only three active neutrinos 
(z/g, u^ and z/^) or more than three light neutrinos. Using A^g to denote the effective 
number of thermally excited neutrino species in the early Universe, some authors have 
recently obtained N^g = 4.34to;88 at the 68% confidence level from an analysis of the 
7-year WMAP data on CMB anisotropies and large-scale structures [9] or A^g = 4.78l};79 
at the 95% confidence level from a similar analysis including the SDSS data on the (DR7) 
halo power spectrum [10] . Moreover, two independent groups have recently found slightly 
higher values of the primordial ^He abundance [TI], implying the presence of additional 
radiation or relativistic particles during the BBN epoch. As argued by Hamann et al. 
in Ref. [12], these results drop an interesting hint that there might exist one or more 
light sterile neutrinos in the Cz/B. Their detailed analysis supports this conjecture and is 
compatible with an interpretation of the LSND [13] and MiniBOONE [H] anomalies in 
terms of three active neutrinos and two sterile neutrinos [15] if the mass scale of sterile 
neutrinos lies in the sub-eV range. Although such an interpretation may have severe 
tension with current disappearance experiments of neutrino oscillations [16] , the situation 
remains so confusing that one might better keep all possibilities open. In particular, 
current cosmological data cannot be used as an argument against the existence of light 
sterile neutrinos [12] and the latter might not necessarily be relevant to current neutrino 
oscillation experiments. This mild standpoint motivates us to examine possible effects of 
light sterile neutrinos in a neutrino capture process to directly detect the Cz/B. 

We find that such light sterile neutrinos can leave a distinct imprint on the electron 
energy spectrum in the capture of relic electron neutrinos by means of radioactive beta- 
decaying nuclei. Considering both the (3 -|- 1) and (3 -|- 2) schemes of neutrino mixing, we 
calculate possible signals of sterile neutrinos relative to active neutrinos, characterized by 
their masses and sensitive to their number densities, in the reaction z/^ + ^H — t- '^He + e~ 



against the corresponding tritium beta decay. Although our numerical results are just for 
the purpose of illustration, we stress that this kind of direct laboratory detection of the 
Cz/B and its sterile component might not be hopeless in the long term. 



2 In the presence of N^ species of light sterile neutrinos, the flavor eigenstates of 



three active neutrinos can be written as 

i 

where a runs over e, fi and r, z/^ is a mass eigenstate of active (for 1 < z < 3) or sterile 
(for 4 < z < 3 + Ng) neutrinos, and V^j stands for an element of the 3 x (3 + N^) 
neutrino mixing matrix. For simplicity, we assume that the light sterile neutrinos under 
consideration do not significantly affect the values of two mass-squared differences and 
three mixing angles of active neutrinos extracted from current experimental data on solar, 
atmospheric, reactor and accelerator neutrino oscillations [IJ . In this assumption we shall 
use Am2i ^ 7.6 x 10"^ eV^ and |Am^J ^ 2.4 x 10"^ eV^ together with 6^2 ~ 34° and 
61^ ~ 10° as typical inputs in our numerical estimates. Depending on the sign of Amf^, 
two mass patterns of three active neutrinos are possible: 



the normal hierarchy rrii < m2 = yTnl + Ar7i2i < ""^3 = V^^i + lAmiil 



the inverted hierarchy m^ < rrii = J ml + |Am3]^| < 7712 = \Jrn\ + |Ar?7,|]^| + Amg 
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In either case the absolute mass scale (m^ or mg) is unknown, but its upper bound is 
expected to be of 0(0. 1) eV as constrained by current cosmological data |9]. Following 
Ref. [I2], we assume the masses of sterile neutrinos {jn^, mg, ■ ■ ■) to lie in the sub-eV 
range. Their mixing with active neutrinos is constrained by current neutrino experiments 
and cosmological data and should be at most of C(O.l) [I5l [16]. To illustrate the effect 
of sterile neutrinos in a neutrino capture process, we shall simply take 9^^ ^ 10° (for 
i > 4) in our numerical estimates. So we have \V^^\ ~ 0.804, \V^2\ ~ 0.542, \V^^\ ~ 0.171 
and \VJ ^ 0.174 in the (3 + 1) scheme; or \VJ ^ 0.792, \VJ ^ 0.534, jl/gl ^ 0.168, 
|Vg4| ~ 0.171 and jV^gl ~ 0.174 in the (3 + 2) scheme. We reiterate that these numerical 
inputs are mainly for the purpose of illustration. 

Let us concentrate on the relic neutrino capture reaction 1^^+ ^H — )■ ^He + e^, since 
its corresponding beta decay ^H — )■ ^He + e~ + u^ is being precisely measured in the 
KATRIN experiment [17]. The capture rate for each neutrino mass eigenstate u^ hidden 
in the z/^ state is given by [HI E] 

AASB^6.5ClKJ'yr-'MCi-S (2) 

in which Q = n^ /{n^,} denotes the ratio of the number density of relic u^ neutrinos 
around the Earth to its average value in the Universe. The standard Big Bang model 
predicts (n ) ~ (% ) ~ 56 cm~^ today for each species of active neutrinos, and this 

i i 

prediction is also expected to hold for each species of light sterile neutrinos if they could 
be completely thermalized in the early Universe q Although possible interactions and 



^The sub-eV sterile neutrinos discussed in Ref. [H] and here should be most likely to stay in full 
thermal equilibrium in the early Universe, provided their mixing with active neutrinos is not strongly 
suppressed [Hj. We are indebted to S. Hannestad for clarifying this point to us. 



oscillations in the early Universe could slightly modify the values of (n ) and (n^^ ), such 
corrections are actually unimportant for our numerical estimates [TJ [19]. In Eq. (2) the 
unit MCi measures one megacurie source of tritium (about 100 g or 2.1 x 10^^ atoms). 
Note that each v^ state yields a monoenergetic electron, whose kinetic energy is given by 
Qa + Ej^ with Qa = Msjj — Ms^^ — m-g ~ 18.6 keV for tritium. Because at least 
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two active neutrinos are already non-relativistic today in the Cz/B (i.e., their masses are 
much larger than T^ ^ 1.945 K or (p^) = 3T^ ^ 5 x 10""^ eV) and the involved sterile 
neutrinos are assumed to have masses of C (0.1) eV, we arrive at T^ ~ Q^ + m^ as a good 
approximation. If the mass of the lightest active neutrino is below (p^) (i.e., it remains 

hot or relativistic today), then its energy can be expressed as E^, f 



(p^)'^ + rnf , which 

is of 0{{pj)) and hence has little effect on the overall electron energy spectrum. Given 
a finite energy resolution in practice, the ideally discrete energy lines of the electrons 
emitted from the reaction z/g + ^H — )■ ^He + e~ must spread and form a continuous 
spectrum. As usual, we consider a Gaussian energy resolution function defined by 
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(3) 



where Tg is the overall kinetic energy of the electrons detected in the experiment. Using 
A to denote the experimental energy resolution (i.e., the full width at half maximum of 
a Gaussian energy resolution for the outgoing electrons [H [7]), we have A = 2v21ii2 a ^ 
2.35482 a . Then the overall neutrino capture rate (i.e., the energy spectrum of the detected 
electrons for the reaction i/g + '^H — )■ ■^He + e~) is given as 



A/'c.B = EA/'SB^lTg, Tl) ^ 6.5 $:0 |V;.pi?(Tg, Tl) yr-^ MCi" 



(4) 



Taking account of the gravitational clustering of relic non-relativistic neutrinos around 
the Earth |20], we expect Cj > 1 in general. For simplicity, we shall first assume Cj = 1 in 
our numerical estimates of the capture rate and then illustrate a possible enhancement of 
the signal due to much larger values of Cj. 

The main background of the neutrino capture process under consideration is the stan- 



dard tritium beta decay H — t- He + e 



■v^ 



In this process the effect induced by nonzero 



neutrino masses can show up near the electron's endpoint energy Qn — min(mj), where 
min(mj) means the lightest neutrino mass among vn^ (for i =1, 2, ■••, S + A^^). The finite 
energy resolution may push the above endpoint towards a higher energy region, and hence 
it is likely to mimic the desired signal of the neutrino capture reaction. Given a finite 
energy resolution described by the Gaussian function in Eq. (3), the energy spectrum of 
the tritium beta decay can be expressed as 
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where E^ 



T' 



m. 



is the electron energy with Tl being its kinetic component, Nj, 



2.1 X 10 denotes the target factor whose value is equal to the number of tritium atoms 



of the target, F(Z,E^) represents the Fermi function, and |A1p ~ 5.55 stands for the 
dimensionless contribution of relevant nuclear matrix elements |T7]. In Eq. (5) the theta 
function 6(Qa — T^ — m-) is adopted to ensure the kinematic requirement. 



With the help of Eqs. (4) and (5), one may numerically calculate the relic 
neutrino capture rate against the corresponding tritium beta decay. Blennow has done 
such an analysis in the scheme of three active neutrinos by examining the dependence 
of both the Ci^B signal and its background on the neutrino mass hierarchy and neutrino 
mixing angles [7]. Here we focus on the sterile component of the Cz/B and its possible 
signals in the neutrino capture reaction i>^+ ^H — )■ ^He + e~. As mentioned in section 2, 
we are mainly interested the sub-eV sterile neutrinos whose mixing with active neutrinos 
is of (9(0.1) or somewhat smaller. To be more explicit, we assume that the masses of 
sterile neutrinos are larger than the absolute mass scale of three active neutrinos. In this 
case we have the following naive expectations: 

• The signal of the sterile component of the Cz/B is on the right-hand side of the 
electron energy spectrum as compared with the signal of the active component of 
the Cz/B. Their separation is measured by their mass differences. 

• The rate of events for the signal of relic sterile neutrinos is crucially dependent on 
the magnitude of their mixing with active neutrinos. Hence a larger value of |V"gj| 
(for i > 4) leads to a higher rate of signal events. 

• Whether a signal can be separated from its background depends on the finite energy 
resolution A in a realistic experiment. In general, A < m-/2 (for z = 1, 2, ■ ■ ■ , 3 + A^^) 
is required to detect the Cz/B via a neutrino capture reaction [HI E]. 

Let us make some quantitative estimates of the Cz/B signals in two schemes of neutrino 
mixing: (a) the (3 + 1) scheme with one sterile neutrino; and (b) the (3 + 2) scheme 
with two sterile neutrinos. In either case the lightest active neutrino mass is typically 
taken to be 0.0 eV, 0.05 eV or 0.1 eV. We fix m^ = 0.3 eV together with \V^^\ ^ 0.804, 
\V^2\ ~ 0.542, iVggI ^ 0.171 and \V^^\ ^ 0.174 in the (3 + 1) scheme; or m^ = 0.2 eV and 
m^ = 0.4 eV together with \V^^\ ^ 0.792, II/2I ~ 0.534, \VJ ^ 0.168, \VJ ^ 0.171 and 
iVggI ~ 0.174 in the (3 + 2) scheme. The gravitational clustering of relic neutrinos around 
the Earth is tentatively omitted and will be illustrated later. Our numerical results are 
presented in Figs. 1 — 3. Some discussions are in order. 

Fig. 1 shows the relic neutrino capture rate as a function of the kinetic energy Tg of 
electrons in the (3 + 1) scheme with A?7i|^ > and m^ = 0.3 eV. The value of the finite 
energy resolution A is taken in such a way that only the signal of the sterile neutrino can 
be seen (left panel) or both the signals of active and sterile neutrinos can be seen (right 
panel). These simple results confirm the naive expectations given above, especially for 
the signal of the sterile component of the Cz/B. We find that it is in principle possible 
to distinguish the sterile neutrino from the background when A is of (9(0.1) eV or much 
smaller. As the lightest neutrino mass m^ increases from to 0.1 eV, the signal curve 
moves towards the higher T^ — Qg region while the background curve moves towards the 
lower T^—Qo region. To make a realistic measurement sensitive to the active component of 
the Cz/B, one needs a sufficiently good energy resolution. The required energy resolution 
depends on the mass hierarchy of three active neutrinos, as one can see from Fig. 1 (right 



panel). Note that the small peak sitting at T^ — Q^ ^ 0.05 eV in the top right corner of 
Fig. 1 arises from the contribution of v^ with m^ ^ 0.05 eV. It is not seeable in the top 
left corner of Fig. 1 simply because the resolution is not good enough. 

An analogous analysis of the relic neutrino capture rate is carried out in the (3 + 1) 
scheme with Am|^ < and m^ = 0.3 eV, and the numerical results are given in Fig. 2. 
We see that this figure is essentially similar to Fig. 1, and the primary difference appears 
in the signal peaks of active neutrinos in their top right corners. Such a difference can 
be understood with the help of Eqs. (3) and (4): the central position of a signal peak is 
approximately determined by T^ — Q^ ~ E,^ and its height is mainly measured by | KjP for 
given values of A and (^. When the lightest neutrino mass m^ is extremely small, the main 
contribution of active neutrinos to the capture rate sits at T^ — Qp^m^'^ m^ ~ 0.05 eV 
as shown in the top right corner of Fig. 2. In contrast, the main signal peaks shows up at 
T^ — Qp G {rrii, mj} when the lightest neutrino mass m^ is extremely small (as illustrated 
in the top right corner of Fig. 1). So it is in practice much easier to detect the active 
component of the Ci^B when three active neutrinos have an inverted mass hierarchy or a 
nearly degenerate mass spectrum. 

We illustrate the relic neutrino capture rate in the (3 + 2) scheme with Arrig^ > 
(left panel) or Arrigj^ < (right panel) in Fig. 3, where m^ = 0.2 eV and m^ = 0.4 eV 
are typically taken. The value of the finite energy resolution A is chosen in such a way 
that only the signals of sterile neutrinos can be seen. If the values of m^ and m^ are very 
close to each other or their difference is much smaller than A, it will be very difficult to 
distinguish between the two signal peaks of sterile neutrinos. When the absolute mass 
scale of three active neutrinos is very close to the smaller mass of two sterile neutrinos, 
there may be a mixture between their signals as shown in the bottom left and bottom 
right corners of Fig. 3. In this case only the heavier sterile neutrino is distinguishable 
from the signal of active neutrinos. 

So far we have neglected possible gravitational clustering of relic neutrinos in our 
numerical estimates. To illustrate this effect, we calculate the relic neutrino capture rate 
in the (3 + 2) scheme with m^ = or m^ = 0, m^ = 0.2 eV and m-g = 0.4 eV. We assume 
a larger gravitational clustering effect for a heavier neutrino around the Earth [20]. In 
this reasonable assumption we typically take Ci = C2 = Cs = 1 (without clustering effects 
for three active neutrinos because their maximal mass is about 0.05 eV in the scenario 
under discussion) and Q = 2(^ = 10 (with mild clustering effects for two sterile neutrinos 
because their masses are 0.2 eV and 0.4 eV, respectively). The value of the finite energy 
resolution A is chosen in such a way that only the signals of sterile neutrinos can be 
observed. As shown in Fig. 4, the signals of two sterile neutrinos are obviously enhanced 
due to C4 > 1 and Q > 1. If the gravitational clustering of non-relativistic neutrinos 
is very significant around the Earth, it will be very helpful for us to detect the Cz/B by 
means of the neutrino capture processes. 



Motivated by a mild standpoint that current cosmological data cannot be used 
as an argument against the existence of light sterile neutrinos and the latter might not 
necessarily be very relevant to current neutrino oscillation experiments, we have examined 
possible effects of sub-eV sterile neutrinos in a neutrino capture reaction to directly detect 
the Ci^B. We find that such light sterile neutrinos can in principle leave a distinct imprint 
on the electron energy spectrum in the capture of relic electron neutrinos by means of 



radioactive beta-decaying nuclei. Considering both the (3 + 1) and (3 + 2) schemes of 
neutrino mixing, we have calculated possible signals of sterile neutrinos relative to active 
neutrinos in the neutrino capture process z/^ + ^H — )■ ^He + e~ against the corresponding 
tritium beta decay. We reiterate that our numerical results are just for the purpose of 
illustration, but they are useful to show a few salient features of the active and sterile 
components of the Cz/B in such neutrino capture processes. 

The major limiting factor of a neutrino capture experiment is its energy resolution. If 
the energy resolution is gradually improved, it should be possible to establish a signal of 
the Cz/B beyond the endpoint of the electron energy spectrum of the beta decay in the 
future. One might only be able to observe one peak for a given energy resolution, either 
because the peaks of different neutrino mass eigenstates merge into a wide one or because 
the peaks of lighter neutrinos are overwhelmed by the background (or for both reasons). 
In this case some information on the sterile component of the Cz/B could be obtained 
from a combined analysis of the available data from the neutrino capture experiment, 
neutrino oscillations and cosmological observations. If multiple peaks could finally be 
seen in the electron energy spectrum with a sufficiently good energy resolution, then they 
would convincingly indicate the presence of light sterile neutrinos in the Cz/B. 

Needless to say, a direct detection of the Cz/B in a realistic experiment is very com- 
plicated and sophisticated. We admit that current experimental techniques are unable 
to lead us to a guaranteed detection of the Cz/B, but more and more efforts have been 
being made towards this ultimate goal. Whether the neutrino capture experiment (e.g., 
v^+ ^H — )■ ^He + e~ under discussion) can win the race or not remains an open question, 
but it is certainly the most promising horse at present. 

One of us (Z.Z.X.) would like to thank S. Hannestad for a very helpful discussion. This 
work is supported in part by the National Natural Science Foundation of China under 
grant No. 10425522 and No. 10875131. 
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Figure 1: The relic neutrino capture rate as a function of the kinetic energy T^ of electrons 
in the (3 + 1) scheme with Am|^ > and m^ = 0.3 eV. The solid and dashed curves 
represent the Cz/B signal and its background, respectively. The value of the finite energy 
resolution A is taken in such a way that only the signal of the sterile neutrino can be seen 
(left panel) or both the signals of active and sterile neutrinos can be seen (right panel). 
The gravitational clustering of relic neutrinos around the Earth has been omitted. 
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Figure 2: The relic neutrino capture rate as a function of the kinetic energy T^ of electrons 
in the (3 + 1) scheme with Am|^ < and m^ = 0.3 eV. The solid and dashed curves 
represent the Cz/B signal and its background, respectively. The value of the finite energy 
resolution A is taken in such a way that only the signal of the sterile neutrino can be seen 
(left panel) or both the signals of active and sterile neutrinos can be seen (right panel). 
The gravitational clustering of relic neutrinos around the Earth has been omitted. 
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Figure 3: The relic neutrino capture rate as a function of the kinetic energy T^ of electrons 
in the (3 + 2) scheme with Am|^ > (left panel) or Arrigj^ < (right panel). In either 
case m^ = 0.2 eV and m^ = 0.4 eV are typically taken. The solid and dashed curves 
represent the Cz/B signal and its background, respectively. The value of the finite energy 
resolution A is taken in such a way that only the signals of sterile neutrinos can be seen. 
The gravitational clustering of relic neutrinos around the Earth has been omitted. 
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Figure 4: The relic neutrino capture rate as a function of the kinetic energy Tg of electrons 
in the (3 + 2) scheme with Am|^ > (left panel) or Amgj^ < (right panel). In either 
case m4 = 0.2 eV and m^ = 0.4 eV are typically taken. The solid and dashed curves 
represent the Cz/B signal and its background, respectively. The value of the finite energy 
resolution A is taken in such a way that only the signals of sterile neutrinos can be seen. 
The gravitational clustering of relic sterile neutrinos around the Earth has been illustrated 
by taking Ci = C2 = Cs = 1 ^^^ C5 = 2C4 = 10 for example. 
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